Vortex behavior near a spin vacancy in 2D XY-magnets 
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The dynamical behavior of anisotropic two dimensional Heisenberg models is still a matter of 
controversy. The existence of a central peak at all temperatures and a rich structure of magnon 
peaks are not yet understood. It seems that the central peaks are related, in some way, to structures 
like vortices. In order to contribute to the discussion of the dynamical behavior of the model we use 
Monte Carlo and spin dynamics simulations as well analytical calculations to study the behavior of 
vortices in the presence of nonmagnetic impurities. Our simulations show that vortices are attracted 
and trapped by the impurities. Using this result we show that if we suppose that vortices are not 
very much disturbed by the presence of the impurities, then they work as an attractive potential to 
the vortices explaining the observed behavior in our simulations. 
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The anisotropic Heisenberg model (AHM) in two di- 
mensions has received a lot of attention recently. Such 
attention is grounded in the fact that a large variety of 
models may be mapped in the AHM. It is an interest- 
ing model because it can support topological excitations, 
like solitons and vortices, which are present in several im- 
portant phenomena. Topological defects are present in 
condensed matter systems such as superconductors, liq- 
uid crystals, superfluids, magnetic materials and several 
others. The knowledge of how such structures behave is 
essential for technological applications and an important 
condition for the understanding of many physical ques- 
tions. Of great importance is the dynamical behavior of 
vortices in ordered structures. Knowing how vortices are 
pinned in superconductors is essential for applications 
such as magnetic levitation, improving magnetic reso- 
nance imaging devised for medical diagnosis and many 
others. Topological defects can be the signature left be- 
hind by the cosmological phase transitions, which occur 
while the universe expands and cools. Beside all of that, 
the study of the dynamics of the AHM model is inter- 
esting by itself. There are several questions not yet re- 
sponded about the model. For example, the origin of 
the central peak in the dynamical structure factor ob- 
served in experiments and simulations is the source of 
several controversial interpretations. Its origin may be 
due to vortex motioni*^. However this point is still 
controversial^. The existence of some kind of impu- 
rity in the system may affect its properties in several 
ways. For example, solitons near a nonmagnetic impu- 
rity in 2D antiferromagnets cause observable effects in 
EPR experiments&t Lattice defects such as impurities 
and dislocations play a crucial role in disrupting order in 



solids. 

The main task of this letter is to consider the dynam- 
ical effects caused by impurities on topological defects. 
We will consider the classical two-dimensional magnetic 
(XY) model described by H = -JE^n^^+M- 
Here J is a coupling constant, S a (a — x, y, z) are the 
components of the classical spin vector S — (S x , S v , S z ) 
and the summation is over nearest-neighbor in a square 
lattice. The spin field can be parametrized in terms of 
spherical angles as S — (cos9cos<f>, cosOsin<f>, sinff), where 
we took \S\ — 1. The continuum version of the XY 
Hamiltonian can be written as 



H=~ I d 2 r 



m 2 (Vm) 2 



(1) 



(l-m 2 )(V</>) 2 + ^m 2 



where m = sinO and a is the lattice constant. The 
spin dynamics is given by the equations of motion, rh — 
8h/5<f>, <j) = —Sh/Sm, where h is the Hamiltonian density 
from Eq.JQ). The introduction of nonmagnetic impuri- 
ties modify the Hamiltonian in the following w&y&&. As 
the model considers only spin interactions up to nearest- 
neighbors, to take into account the presence of a non- 
magnetic impurity we introduce a circular hole with ra- 
dius of the order of the lattice spacing, inside of which 
the Hamiltonian density vanishes. The Hamiltonian, Eq. 
||TJ, becomes^ 



m 2 (Vm) 2 



(1 -TO 2 )(V0) 2 
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(2) 



U(r) 
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with U(r) — 1 if | r — r |> a and C/(r) = if | r — r |< a, 
for a nonmagnetic impurity at r . Early works&S have 
found a repulsive force between the vortex and the non- 
magnetic vacancy. The principal ingredient in the calcu- 
lations was to suppose that the spin vacancy defforms 
considerably the vortex structure. However, although 
the assumption of vortex deformation near the impurity 
seems reasonable, it needs a convincing demonstration. 
We discuss this question by means of numerical simula- 
tions, which will be used not only to check the validity 
of the approximations done, but also as a guide to give 
us a clue on which assumptions can be made in a safe 
way. We have performed spin dynamical simulations on 
a L = 20 square lattice to observe the behavior of a single 
vortex initially located at the center of the system, in the 
presence of a spin vacancy located two sites away from 
the center of the system. In order to introduce a sin- 
gle vortex into the system, we have imposed diagonally 
antiperiodic boundary conditions^ 



Sl+i, v 
S x ,l+\ 



— —Si.L-y+l ! So y — - 



~Sl-x+1,L, 



(3) 



for all 1 < x, y < L. The discrete equation of motion for 
each spin isii 



dS m 
dx 



ef 



where 



H ef = -jJ2(SZe x + SV n e y ) 



(4) 



(5) 



and e x and e y are unit vectors in the x and y directions, 
respectively. The equations of motion were integrated nu- 
merically forward in time using a vectorized fourth-order 
Runge-Kutta scheme with a time step of 0.04J -1 . The 
initial structure is shown in Fig^ In contrast to the an- 



toward the spinless site, indicating an effective attractive 
potential of interaction between the vortex and the va- 
cancy, as shown in Figs [3 and |21 In Fig|21we show the 
configuration after 70 time steps and Fig|3|after 150 time 
steps. We notice that after 150 time steps the position of 
the vortex center reaches the equilibrium with its center 
at the spinless site. The key point here is the property of 
non-deformation of the vortex configuration. The attrac- 
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FIG. 2: Configuration after 70 time steps. 
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FIG. 3: Configuration after 150 time steps. 
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tive potential could also be obtained from our theoretical 
model Hamiltonian (J2J) by assuming that the vortex con- 
figuration is not deformed due to the presence of an im- 
purity. Substituting both the non-deformed, static and 
planar vortex solution 4>q = arctan(y / x) located at the 
origin and a nonmagnetic impurity placed at the site fo 
in Hamiltonian we obtain 



(G) 



iy^fd 2 x~ I (V^ ) 2 cL4(fo) 

A(r ) 



FIG. 1: A typical initial configuration. The square represents 
the spinless site. 

alytical results presented in Ref£, the simulation results 
show that the vortex structure does not change appre- 
ciably near the spin vacancy. The vortex center moves 



where A(fo) is the area of the hole with radius a around 
the point fo that represents the spinless site. The inte- 
grals are easily calculated leading to 



E I =E v + ^ln(l-^ 



(7) 
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where Ei^Tq) and E v = nJln(L/ao) are the vortex en- 
ergies in the presence and absence of the nonmagnetic 
impurity respectively for a system of size L. The con- 
stant ao — 0.24c 12 present in E v> is some suitable short- 
distance cutoff to avoid spurious divergences due to the 
fact that the vortex center is a singularity in the contin- 
uum limit. Its appropriated value for the square lattice 
was determined numerically in Ref^. Of course, Eq.Q 
is not valid in the limit ro — > 0, where the continuum ap- 
proximation breaks down. In this limit there will have an 
intersection between the circular hole of the impurity and 
the structure of the vortex core. A point of the circumfer- 
ence of the impurity circle would meet the center of the 
vortex core lowering the vortex energy as ro decreases. 
Then, from Eq.JJJ) the effective potential experienced by 
the two defects can be written as 

V eff (r Q ) = E I (r )-E u (8) 

ttJ, / a 2 \ 
= —m yl- -j j , for r > a, 

which is attractive. For large distances of separation 
(ro >> a), the Eq.© can be approximated by V e ff(r ) rs 
— (ttJo 2 /2)(1/tq). To know the behavior of the potential 
in the region inside the core (ro < a), we first put the hole 
exactly at the vortex center, obtaining Ej — TvJln(L/a), 
which leads to V eff (Q) = £/(0) - E v = itJln{0M) = 
—4. 48 J. Then, using this result, we modify Eq. JHJ) to get 
a complete function V e //(ro), introducing a small con- 
stant b as follows 

Wn.) = ^i»(i-^), (9) 

with the condition 

VeffiP) = ^-In (l - ^) = TJZn(0.24). (10) 

This artifact leads to an expression valid in all space that 
reproduces correctly the asymptotic limits for b — 1.03a. 
The potential (Eq. [§J| indicates that vortex can be 
trapped by a nonmagnetic impurity. The attractive po- 
tential experienced by the vortex can be written as 

where X is the distance between vortex center and the 
spinlcss site. It has a minimum at X — and decreases 
rapidly for X > 2a (See FigH). 

In order to observe the lowest energy configurations at 
low temperatures in a classical two-dimensional system 
described by the XY model, with a nonmagnetic impu- 
rity, we performed Monte Carlo (MC) simulations. The 
simulations were done on a L — 20 square lattices at tem- 
perature T = 0.1 in units of J/ks (kg is the Boltzmann 
constant), using a standart Metropolis algorithm 13 with 
diagonally antiperiodic boundary conditions (Eq[3|) and 
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FIG. 4: Attractive potencial as a function of X/a 



random initial configurations. Using these boundary con- 
ditions, at low temperature, the equilibrium configura- 
tions can have a single vortex or a single anti- vortex. We 
have observed that after 10 5 MC steps the equilibrium 
configurations was reached. Fig|5] shows a single anti- 
vortex equilibrium configuration after 2 x 10 5 MC steps 
with the nonmagnetic impurity located at site (12,10). 
Fig. [S] shows a single vortex equilibrium configuration 
after 2 x 10 5 MC steps with the nonmagnetic impurity 
located at site (13,12). We have observed that, in both 
cases, in the lowest energy configurations, the anti- vortex 
and vortex centers are located at the nonmagnetic impu- 
rity, which is in agreement with Ea. lll|) . 
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FIG. 5: Anti- vortex configuration after 2 x 10 5 MC steps. The 
impurity is located at site (12,10). 

In summary, using Monte Carlo and molecular dynam- 
ics simulations, we observed that a vortex structure put 
close to a non magnetic site in a square lattice does not 
change appreciably its geometry. As the system evolutcs 
in time the vortex center moves toward the spinless site. 
By using a continuum approach we have modeled the in- 
teraction as an attractive potencial which has a minimum 
for the vortex center located at spinless site, as observed 
in MC simulations. The study of the dynamics of diluted 
models can help us to understand the origin of the cen- 
tral peak found in early simulations. The central peak is 
believed to be due to vortex motion. If that is the impor- 
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FIG. 6: Vortex configuration after 2 x 10 5 MC steps, 
impurity is located at site (13,12). 



The 



tant contribution, the presence of non magnetic sites will 
diminish the high of the central peak due to a trapping 
of vortices at that non magnetic sites. 
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